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ABSTRACT 

Context. Red super-giant (RSG) stars exhibit significant mass loss through a slow and dense wind. They are often considered to be the 
more massive counter-parts of Asymptotic Giant Branch (AGB) stars. While the AGB mass-loss is linked to their strong pulsations, 
the RSG are often only weakly variable. This raises the question whether their wind-driving mechanism and the dust composition in 
the wind are the same. 

Aims. To study the conditions at the base of the wind, by determining the dust composition in a sample of RSG. The dust composition 
is thought to be sensitive to the density, temperature and acceleration at the base of the wind. To compare the derived composition 
with the composition found in AGB star winds. 

Methods. We compile a sample of 27 RSG infrared spectra (ISO-SWS) and supplement these with photometric measurements to 
obtain the full spectral energy distribution (SED). These data are modelled using a dust radiative transfer code, taking into account 
the optical properties of relevant candidate materials. The results are scrutinised for correlations in terms of mass-loss rate, density at 
the inner edge of the dust shell and stellar parameters. 

Results. We find (1) strong correlations between dust composition, mass-loss rate and stellar luminosity, roughly in agreement with 
the theoretical dust condensation sequence, (2) the need for a continuous (near-)IR dust opacity and tentatively propose amorphous 
carbon, and (3) significant differences with AGB star winds: presence of PAHs, absence of 'the' 13 /jm band, and a lack of strong 
water bands. 

Conclusions. Dust condensation in RSG is found to experience a similar freeze-out process as in AGB stars. Together with the 
positive effect of the stellar luminosity on the mass-loss rate, this suggests that radiation pressure on dust grains is an important 
ingredient in the driving mechanism. Still, differences with AGB stars are manifold and thus the winds of RSG deserve separate 
studies. 

Key words, stars: circumstellar matter, stars:supergiants, stars:winds, outflows, stars:atmospheres, infrared: stars 



o 

O ! 1. Introduction dHoefner et al][l998t and further updates). Instead, convection, 

L * ■ „ . . . , „ , , . _ „ , „. , . rotation and/or chromospheric activity may turn out to play an 

> ■ As stars of initial mass 8 < M ; lt < 40 M evolve off the main se- . ^ mle For e j JosseHn & ple j ^ ) find a coj ._ 

£ ; quence and reach the core helium-burning phase they turn into reMon between mass . loss rate and me str m of the hoto . 

r> , red super-giant (RSG) stars, the largest and most luminous of all SDner j c turbulence 

■ stars (for a recent review on massive star observation and evo- 

' lution, see lMassevll2003l) . Through their stellar wind and even- While the mechanism responsible for the initial levitation 

tually when they end their life through a supernova explosion, of the atmosphere may therefore be different in RSG as com- 

they enrich the ISM with heavy elements and large amo unts of pared to AGB stars, both types of stars are often used together 

kinetic energy, possibly trigg ering new star formation (Gehri in studies of the du st condensation sequence and wind driv - 

Il989t iHerbst & Assousalll977l) .The dusty stellar winds of RSG ing mechanis m (e.g. iMatsuura et al.ll2005t iDijkstra et aDl2005h . 

not only influence their evolution, but also complicate the deter- Nevertheless, Spec k et al.l(l2000l) find differences in dust compo- 

mination of general properties such as the effective temperature sition between RSG and AGB stars, with the former often show- 

and surface gravity, e.g. throug h reddening which a ppears differ- ing Ca-Al-rich silicates instead of magnesium silicates, 
ent from that through the ISM dMassev et al.ll2005l) . 

Most RSG are irregular variables with smaller amplitudes Assuming a dust_condensation sequence such as that pre- 

and higher effective temperatures than AGB stars. Moreover, sented by | Tielens | dl990|) and placing it in an expanding stellar 

they often show significant chromospheric activity. These prop- Wlnd ' one ex P ects a P rocess called "freeze-out", where the se- 

erties make them intrinsically very different from their lower- 1 uence 18 not completed but stops at an intermediate product 

mass counterparts in which strong pulsations and cool molecu- when the Wlnd denslt y dro P s below the denslt y required for the 

lar layers are believed to be crucial ingredients for the mass loss next step This freeze-out is now tentatively observed in AGB 

stars (e.g. lHeras & Honyl2005l) . If the dust condensation and the 

Send offprint requests to: T. Verhoelst, e-mail: wind dynamics are similar for the RSG, one may also expect to 

tijl.verhoelst@ster.kuleuven.be observe incomplete dust condensation. We test this hypothesis 

* Postdoctoral Fellow of the Fund for Scientific Research, Flanders through the analysis of all ISO-SWS spectra of RSG. 
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In Sect.|2] we present the construction of the sample and its 
main characteristics. Sect.[3]deals with the reduction of the ISO- 
SWS spectra. In the following section, we discuss the details 
of our modelling strategy, including the photosphere models, 
dust components, radiative transfer method and fitting strategy. 
In Sect. 15. II we analyse an "averaged" RSG dust spectrum to 
check for any major short comings in our modelling approach. 
The results on all individual sample stars are then summarised 
in Sect. l5.2l and trends and correlations between the derived pa- 
rameters are presented. The discrepancies between models and 
observed spectra are studied in Sect. [6] We end with a discussion 
and the conclusions (Sect.|8]l. 

2. The sample 

The availability of an ISO-SWS spectrum is an essential ingre- 
dient for our research aims and the first step is therefore the se- 
lection of all stars with spectral types M and K and luminosity 
classes I and II fr om the list o f usable 2.4^-5 /mi ISO-SWS spec- 
tra presented by ISloan et alj (120031) . From this list of 41 stars, 
we remove the RV Tauri star RSct as it is a Post-AGB object. 
The spectra of XX Per, IRC +60370 and WX Cas show technical 
anomalies and are also removed from our list. D istances are de- 
rived from the hipparcos parallax measurement (Perrvma n et alj 
11997b or from the assumed distances to the clusters containing 
our program stars dHumphrevslll978l) . For IRC+40427, no reli- 
able distance estimate is available. We assume it to be a typical 
2 kpc, which may in fact be to nearby, as the observed redden- 
ing is much stronger than the combined effect of that distance 
and its circumstellar shell (see Sect l5.2b . As the strong redden- 
ing may also be a consequence of a non-spherical distribution of 
the circumstellar dust, we refrain from deriving a distance from 
the observed reddening. 

Effective temperatures and bolometric corrections are esti- 
mated from the spe ctral type based on the relations found by 
iLeyesque et al.1 d2005l). Compa rison with the evolutionary tracks 
of iLeieune & S chaerer (2001) reveals 24 bona fide red super- 
giants, and 3 strong candidate^]. VY CMa is an RSG experienc- 
ing extreme mass loss resulting in an optically thick circumstel- 
lar environment (CSE). The detailed analysis required to derive 
reliable dust composition parameters for VY CMa i s beyo nd the 
scope of this paper, and we refer to Har wit et al.1 d200ll) for a 
study of its ISO-SWS spectrum. 

We use these 26 stars for our study. Their designations and 
main characteri stics are presen t ed in TableQ] The o verlap with 
the samples of ISvlvester et alj d 19941) . ISpeck et alj d2000t) and 
Levesque et al. (2005) is 8,1 1 and 14 stars, respectively. 

iKraemer et al.1 d2002) present a classification of, amongst 
others, our sample stars on the basis of their photospheric and 
dust characteristics. Most stars are in the 2.SEa, 2.SEb or 2.SEc 
classes, meaning they show strong dust features from silicates 
and alumina, and molecular bands in the photospheric part of 
the spectrum. 

3. The observations 

The Infrare d Space Observato ry Short Wavelength Spectrometer 
(ISO-SWS iKessler et al.l 1 1 996h observations were done in the 
AOT1 observing mode, resulting in a low -resolution full grat- 
ing scan, except for S Per, which was observed at full instru- 



mental resolution in the AOT06 mode. The observations were 
then processed using the SW S interactive analysis product, IA 
(see Ide Graauw et alj 1 19961) using calibration files and proce- 
dures equivalent to pipeline version 10.1. Further data process- 
ing consisted of extensive bad data removal and rebinning on 
a fixed resolution (/l/A/l=200) wavelength grid. The bad data- 
removal uses the large amount of redundancy in the available 
spectral scans to identify cosmic hits, sudden changes in dark- 
current and other artifacts that cause the signal of a single or 
multiple detectors to diverge from the mean of the other spec- 
tral scans. In order to combine the different subband into one 
continuous spectrum from 2 to 45 /mi we have applied scal- 
ing factors or offsets. In general the match between the different 
subbands is good and the applied scaling/offsets are small com- 
pared to the flux calibration uncertainties with a few exceptions: 
Alpha Her and NR Vul show much larger differences between 
the different subbands than can be expected on the basis of the 
flux-calibration uncertainties alone. This is most likely due to 
mispointed observations. For these two observations we use the 
IRAS point-source fluxes as a guideline for splicing the spec- 
trum. The RMS noise increases with increasing wavelength, in 
particular at A > 27 /mi. The flux levels of our sample stars drop 
with wavelength, in particular if there is little on no dust excess. 
This causes some of the longest wavelength parts of the spec- 
tra to be noise dominated. See as the worst example H 14580 in 
Fig. [2] The noise-dominated parts are not used in the analysis. 

For a proper estimation of the photospheric irradiation 
and underlying continuum, we prefer to construct a spec- 
tral energy distribution (SED) covering wavelengths from 
the UV to the far-IR. Optical pho t ometry can be found i n 
Johnson et al. ( 1966a), Johnson et al. ( 1966b), Mendoza ( 196"/ 



Led d 19701). ICousins & Lagerweiil d 197 lb. iHumphrevs & Nev 
dl974). IWawrukiewicz & Led dl974l) . [Nicoletl dl978l) and 



Wawrukie wicz & Lee 
iKharchenkol d2001l) . Near-IR photome try is available for a ll our 
sources from the 2MASS catalogue (Skrut skie et al.1 120061) and 
the far-IR side of the SED can be covered with the IRAS ob- 
servations which are ava ilable for all but seven of our sources 
dNeugebauer et alJll984b . 

Using interstellar e xtinction estimates from either the model 
bv lArenou et alj d 19921) or from the observed extinction towards 
the early-type stars in the corresponding clusters, we deredden 
the photometric and spe ctroscopic data with the law presented 
by [Cardelli et al. ( 1989), e xtended toward longer wavelengths 
with the local ISM curve of Chiar & Tielensl d2006l) . and assum- 
ing R v = 3.6 dMassev et al.H2005b . 



1 Stars removed from the list because they are either on the AGB 
or even on the RGB: TCet, WXCas, BD+59 594, i Aur, HD90586, 
RT Car, BC Cyg, CIT 1 1 and IRC+60370, p Cap, S Pav and a UMa 



4. Modelling 

The focus of this study is to derive acurate dust masses and its 
composition. Therefore, several aspects of the complex RSG at- 
mospheres are not taken into account and the following simpli- 
fications are made constructing the radiative transfer models: 
(1) all calculations are done in ID-spherical geometry, which 
means we neglect inhomogeneities, (2) we do not take the 
chromosphere into account, (3) we do not include any extra- 
photospheric molecular layers and (4) we assume the mass loss 
to be constant. 



4.1. Photospheres 

The stell ar photospheres are rep resented by marcs (edition 1998) 
models dGustafsson et al.|[2008l) as these are specifically devel- 
oped for cool stars, with an emphasis on molecular opacities 
and the effects of extended atmospheres which do not allow a 
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Table 1. Object designations and general characteristics. Star with their names printed in italic have luminosities which also allow an 
AGB star iden tification. The b olometric magnitudes were computed from the dereddened K band magnitudes, using the bolometric 
corrections of lLevesque et alj d2005l) . 
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a Combined General C atalogue of Variable Stars (Kh olopov et al.lfT998h . 

b iKraemer et all (2002) 

Distances below 1 kpc come from parallax measurements with hipparcos (Perryman et al. 1997). Distances equal to 2 kpc are actually un- 
known and distances greater than 1 kpc are cluster distances from lHumphrevsl dl978l) . yu Cep is an exception: we use the distance of Hu mphreys! 
i ll 9781) instead of the poor parallax measurement. 

d ISM extinction estimates for stars within 1 kpc come from the galactic model of Are nou et al 1l fl992h . All other stars and /j Cep have estimates 
based on the extinction toward early-type stars in the corresponding clusters. 



plane-parallel approximation. They were computed in spherical 
geometry with solar abundances, a surface gravit}0 of logg = 
0.5, a microturbulent velocity v tur b = 2kms~ I and tempera- 
tures ranging from 3000 to 4500 K. The effective temperatures 
were derived from the spectral type following the relations from 
iLeve sque etal.1 (120051) . 

4.2. Dust radiative transfer 

The dust shell is modelled u sing the proprietary spherical ra- 
diativ e transfer code modust dBouwman et alj|2000t [Bouwrnan 
1200 lb . Under the constraint of radiative equilibrium, this code 
solves the monochromatic radiative transfer equation from 
UV/optical to millimetre wave l engths using a Feautrier type so- 
lution method dFeautrierll 1 964t iMihalasll 1 9781) . The code allows 
to have several different dust components of various grain sizes 
and shapes, each with its own temperature distribution. 



2 Lower gravity models failed to converge at the lowest temperatures 
which led us to adopt this value uniformly. The stellar mass of the mod- 
els in this grid is only 1 M , but at the resolution of our set of observa- 
tion, the difference with higher-mass models would not be noticeable. 



4.3. Dust composition 

We searched the literature for r elevant dust species to include 
in our study. ISpeck et"aT] d2000t) find in the RSG in their sam- 
ple melilite, alumina and olivine^ , but they do not need a 
13 yum-carrier candidate such as Spinel as they do not detect this 
feature in thei r RSG spectra and neither do we (see Sect.|6]l. 
iHeras & Honvl d2005l) find i n thei r sample of AGB stars those 
species used by Speck et al. (2000) and also Mgo.1Feo.9O, which 
has a feature at 1 9.5 /um , outside the w avelength band observed 
by ISpeck et all J2000h . ICarnl 0002) uses olivines, alumina, 
spinel and Mgo jFeo.ciO. For a discussion as to the need for 
m elilite instead of a n iron-magnesium silicate in RSG, we refer 
to lSpeck et alj d2000l) . who show the existence of a category of 
RSG 10 micron spectra with a peak position significantly to the 
red of that of olivines. Moreover, they show an excellent agree- 
ment with a predicted melilite emission spectrum. 

From a theoretical point of view, a dust formation scenario 
usually starts with TiC>2 seeds, on which first simple oxides such 



3 Remark that the term Olivine in principle refers to the crystalline 
form. Although the correct phrasing should be "amorphous silicate with 
an olivine stoichiometry", we refer to this material as Olivines, as is 
done in similar studies. 
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as AI2O3 and MgFeO grow. If the density/temperature condi- 
tions are right, also melilite, and olivines and pyroxenes with 
iron content can condense onto the grains. The full condensation 
scheme consists of quite a few more possible dust species (e.g. 
Anorthite and Plagioclase), but no spectral sign of these dust 
types has so far been seen in RSG. 

The dust species used in this study and their relevant char- 
acteristics are presented in Table[2] The optical constants are ob- 
tained from the AIU Jena databasfl 

The RSG in this sample appear to be pure outflow sources, 
i.e. nothing indicates that they may be binary systems with sta- 
ble circumbinary dust disks: as shown in Fig. [2] there is no 
strong discrepancy between reddening and IR excess, nor do 
we see any sign of processed, crystalline g rains such as those 
seen around some Post-AGB objects (e.g. Gielenetal. 2007) 
or in protoplanetary disks (e.g. Ivan Boekel et alj [2005b . This 
implies that grain growth beyond 0.1 nm is unli kely (see e.g. 
iGail & Sedlma^ll999tlWoitke & Niccolinill2005l) . The absorp- 
tion cross sections can therefore be calculated in the Rayleigh 
limit, and an underlying grain size distribution is only required 
to derive the dust m asses. We use the size distribution derived 
bv lMathis et al.l ( 1 1977b and compute the cross section in the ap- 
proximation of a Contin uous Distribution of Ellipsoids (CDE, 
iBohren & Huffmadfl983l) . as we found this to provide a better 
agreement with the observations than the Mie approximation. 
More sophisticated treatments of grain sizes and shapes exist 
now, such as a Distribution of Hollow Sphe res (DHS), but w ithin 
the Rayleigh limit, CDE works very well dMin et al.ll2005l). We 
must p oint out that according to very recent work by iHofnerl 
(2008), relaxing the restriction to the small particle limit in the 
modelling of AGB winds leads to a possible grain growth up to 
1 yizm, which may also help in the radiative driving of the wind. 

An overview of the spectral features of these different dust 
species is shown in Fig.Q] An in-depth discussion of this figure 
is presented in Sect. 15. II 

4.4. The strategy 

The unknown angular scale of the central star is determined 
by scaling the model to the ISO-SWS spectrum between 3.5 
and 4.5 yum as these wavelengths are relatively free of dust and 
molecular emission/absorption. This method allows us to bypass 
the uncertain absolute luminosity and distance. The geometrical 
properties of the dust shell are proportional to the stellar radius 
and thus their angular scale is now also fixed. Under the assump- 
tion of a continuous wind, the outer radius of the dust shell can 
not be determined from the current modelling as the outer re- 
gions of the wind are too cold to contribute to the Mid-IR spec- 
trum. We therefore fix its value at 5000 R*. 

The mass-loss rate does depend on the assumed distance and 
luminosity. The initial model grid (central star + dust shell) is 
calculated for a typical luminosity of 10 5 L . This model mass- 
loss rate M gr id is converted to an actual M using the distances 
presented in Sect.|2]and the angular scale derived from the fitting 
(see also Sect. l5.21 >. The outflow velocity and dust-to-gas ratio 
are assumed to be 10km s _1 and 0.01 respectively, both typical 
values for RSG. Note that these are uncertain to a factor of at 
least 2. 

The parameters to be determined from the individual fitting 
of the spectra are the inner radius, R; n , the mass-loss rate, M, and 
the relative abundances of the different dust species. 



http://www.astro.uni-jena.de/Laboratory/Database/databases.html 
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Fig. 1. The average of all dereddened observed dust spectra, each 
normalized with its photospheric flux at 3.5 fim. For comparison, 
the pure photospheric spectrum of a UMa is shown as a black 
dotted line. The spectral features of the commonly used dust 
species are also shown (blue) in the way they appear in models 
with CSEs containing only that dust type. None of these dust 
species can explain the excess already observed before 8 fim. 
The red line is a model containing also some Fe and/or C. 



We refrain from doing a formal minimisation while deter- 
mining the optimum dust-shell parameters. Instead we deter- 
mine the best-fit model using visual inspection, following a well 
defined strategy. The reasons for doing so are two-fold. 1) It is 
our aim to get, at the same time, a good fit to the SED and the 
detailed IR-spectrum. This is a requirement as the circumstel- 
lar reddening puts a strong constraint on the inner shell radius, 
which in turn determines the relative contribution of the individ- 
ual dust species to the observed features through their different 
temperatures. It is virtually impossible to devise a good weight- 
ing scheme to give equal importance to the ~10 photometric 
points and the thousands of spectral points in the S WS spectra. In 
particular, since the uncertainties on the data vary considerably 
throughout the sample, this would imply adapting the weight- 
ing on a per source basis. The only way to judge whether such 
a complicated weighting scheme works well would be by visual 
inspection. 2) There are some spectral regions that may be poorly 
reproduced even with a, to our judgement, good model. This is 
due to the possible presence of components in the real spectra 
which are lacking in the model, like extended molecular layers 
(also known as MOLspheres) or PAHs. Again, such regions and 
there importance can only be identified by visual inspection. 

The core of our strategy is an iteration of the following se- 
quence: 

- We estimate the dust mass-loss rate from the full SED and 
the observed Lir 



5 Lir is calculated as the integration over A > 2 /jm of the final 
model spectrum minus the photospheric model and L* is calculated 
by integrating the photospheric model over all wavelengths. 
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Table 2. Dust species used in this stud y with their commo nly-used name, their size, shape and the source of the optical constants. 
The size distribution is that derived by Mat his et alj d 1977b and the cross section are computed in the CDE (Continuous Distribution 
of Ellipsoids) approximation. 



Name 


Composition 


Lattice structure 


Size (N( 


a)~ or 3 - 5 ) 


Shape 


Reference 


Alumina 


A1 2 3 


amorphous 


0.01 


- 1/jm 


CDE 


Begemann et al. (1997); Koike et al. (1995) 


Melilite 


Ca 2 Al 2 Si0 7 


amorphous 


0.01 


- l/im 


CDE 


Mutschke et al. (1998); Jaeaer et al. (1994) 


Olivine 


Mg .gFei. 2 SiO4 


amorphous 


0.01 


- l/im 


CDE 


Dorschner et al. (1995) 


MgFeO 


Mgo.1Feo.9O 


amorphous 


0.01 


- 1/un 


CDE 


Henning et al. (1995) 


Metallic iron 


Fe 


crystalline 


0.01 


- 1/jm 


CDE 


Henning & Stognienko (1996) 


Carbon 


C 


amorphous 


0.01 


- 1/jm 


CDE 


Preibisch et al. (1993) 



- The relative abundances of melilite, olivine and alumina are 
adjusted to match the peak position of the 10 /jm feature 

- and the same is done for the 18 /mi feature peak position 
using olivine, melilite and Mgo.1Feo.9O. 

- When the peak positions of the features are well reproduced 
but not yet their strength, we vary the inner radius to adjust 
their relative strenghts and the mass-loss rate to alter the ab- 
solute emission levels. Meanwhile, we check that our model 
predicts the correct amount of circumstellar reddening. 

5. Results 

5.1. A continuum opacity source 

In Fig. Q] we present an "averaged" RSG ISO-SWS spectrum, 
where each individual spectrum is dereddened with the appropri- 
ate value (last column in TableQ]) and normalized to the average 
flux between 3.5 and 4 /mi. We exclude RSG without signifi- 
cant mass loss. Also shown are models with envelopes contain- 
ing only a single dust type. 

Striking in this figure is that the observed excess appears to 
start already at 5 /mi, while the modelled dust features are all 
located beyond 8 //m. A related problem could be the slope of 
the photospheric(?) part of the modelled spectrum at 2.5-3 /mi. 

The excess around 6 //m and the difference in slope around 
3 /mi have recently been attributed to the presence of a 
MOLsphere, i.e. an extr a-photospheric layer of molecular ma- 
terial (e.g. iTsuiil I2000), in which especially H2O would be 
the main source of excess emission and absorption. Individual 
lines of water in these MOLspheres have been detected (e.g. 
iJennings & Sadalll998l) and the she lls have been resolved by 
near and mid-IR interferometry (e.g. IPerrin et al.|l2004t lOhnakal 
12004 IPerrin et all 120051) . I Verhoelst et all (120061) argue that a 
purely molecular MOLsphere cannot explain the wavelength de- 
pendence of the opacity across the entire IR range in the case of 
RSG a Ori and they suggest alumina as additional mid-IR opac- 
ity sourc e. However , this d oes not yet explain the observation 
made by iRvde et all d2006l) in the spectra of // Cep, that water 
lines appear in absorption at 12 /an, while they are assumed to 
generate an excess at 6 /itrQ. Moreover, the excess around 6 /mi 
in the individual stars does not show any water lines in propor- 
tion to the excess w.r.t. the photospheric continuum. 

We propose here that, although extra-photospheric molecu- 
lar material is clearly present, the peculiarities in the shape of 
the overall SED discussed here are due to a source of contin- 
uous opacity with a fairly cool temperature. Free-free emission 
from a chromosphere or ionised wind can be present in RSG 

6 The MOLspheres are found to be colder than the photospheric 
background, and generate excess emission only through their larger 
emitting surface. 



(e.g. lHarper et al. 2001), but since its source function has a tem- 
perature above that of the photosphere, it can not explain the ad- 
ditional extinction towards the near-IFQ. The color temperature 
of the excess therefor points to dust rather than free-free emis- 
sion. We find that the inclusion of either metallic Fe, amorphous 
C or micron-sized grains in the dust shell can explain the ex- 
cess emission at 6 /mi, the slope issue at 3 /mi and the additional 
extinction at near-IR wavelengths (see Sect. l5.21 i. The remaining 
residuals are molecular absorption, as expected from the pres- 
ence of a cool MOLsphere. 

Whether this missing dust species is either metallic iron, 
micron-sized grains or amorphous carbon, can not be fully an- 
swered with the current analysis. The presence o f metallic Fe 
in the CSE of VY CMa has been postulated by lHarwit et all 
(120011) and it is als o used to explain the SEDs of OH/IR stars 
(Ke mper et al.ll200"2b . Moreover, metallic Fe is a product of the 
O-rich condensation sequence. Nevertheless, amorphous carbon 
is an equally interesting possibility, as it would be roughly 20 
times as eff ective in accelerating the wind as Fe (see Table 1 in 
Woitke 2006). The occurence of pure carbon in an oxygen-rich 
CSE is plausible dHofner & Andersen] |2007|) . even more so in 
RSG where the chromospheric radiation field can be responsible 
for the dissociation of CO. We find that, for the amount of pure 
carbon required to explain the RSG spectra, only of the order 
of 1 % of the CO must be dissociated. This is well below what 
is found po ssible in theoreti cal models of M stars with chro- 
mospheres dBeck et al.| [l992). More evidence for the presence 
of carbon not bound in CO is the detection of PAH emission 
in four of our sources (Sect. 16. 2731 . Papers discussing the possi- 
ble presence of amorphous carbon in oxygen-rich outflows often 
dismiss the spectroscopic detection of its presence as impossible 
due to the lack of spectral features. Interestingly, we find here 
that, when studying the full SED, it is possible to detect con- 
tinuous opacity sources. Continuous opacity by micron-sized O- 
rich grains, playing a crucial rol e in the driving of the wind, has 
been proposed for AGB stars by Hofnet] d2008l) . but it is unclear 
whether this grain growth can also occur in RSG. 

Regardless of the origin of this continuum emission, it must 
be included in our modelling to derive reliable mass fractions for 
the other dust species. We have chosen to include amorphous C, 
in the knowledge that Fe would have a similar effect. 

5.2. Models for all sample stars 

Figures|2] and [3] present the observed SEDs and spectra respec- 
tively and the best-fit models, sorted following increasing dust 



7 We know that the near-IR is suffering additional extinction (as opp- 
posed to less excess emission) from the fit of our photosphere models 
to the unreddened optical fluxes 
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1 10 1 10 

Wavelength [jum] Wavelength [yu.m] 



Fig. 2. The SEDs of all sample stars with their best-fit models (total spectrum in red dashes and non-reddened photospheric spectrum 
in bluedots) The photometric points are either the 7 Geneva filters or B and V in the Johnson system, J H and Ks from 2mass and 
112 and 125 from the IRAS mission. The dotted vertical line indicates 9.7 jum, the wavelength where "classical" amorphous silicates 




Fig. 3. The ISO-SWS spectra of all sample stars with their best-fit models (symbols and colours as in Fig. [2]). 
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luminosity or following decreasing photospheric temperature in 
case no dust emission is present. In general the agreement is ex- 
cellent. IRC +40427 appears much more reddened than our esti- 
mate of Ay = 1.79, which suggests it to be located much further 
away than our default value of 2 kpc. 

The results on all 27 sample stars are summarized in Table[3] 
For 6 objects, the mass-loss rate is too low (M < 10 _9 M G yr _1 ) 
to allow for a reliable determination of the dust composition. 
The distance and luminosity corrected mass-loss rate presented 
in column 9 is computed as 

M = M grid A(d/lkpc) 2 

where M gr jd is the mass-loss rate for an assumed luminosity of 
10 5 L o (the default value in our model grid), A is the scaling 
factor used to match the observed flux around 4 /urn and d is the 
distance. 



6. Analysis of the residuals 

A general problem, especially for the higher mass-loss stars, is 
the over-prediction of the near-IR photometry. This corresponds 
to the need for an extinction law with a higher Ry. The dust al- 
ready includes a component with grey extinction properties (see 
Sect. |5.11 l, but in our models it is located too far from the cen- 
tral star to provide sufficient extinction. Our model assumes the 
same inner radius for all dust species, but the discrepancy ob- 
served here may indicate that the continuous opacity source is 
located even closer to the base of the wind, right above the pho- 
tosphere. The d etection of alumina 0.5 above the pho tosphere 
of Betelgeuse dVerhoelst et al.1 120061: iPerrin et all 120071) shows 
that dust condensation so close to the photosphere is possible. 

6.1. Molecular bands 

The ISO spectra up to 8 pm are dominated by photospheric and 
possibly continuum dust emission. The observed spectral fea- 
tures should therefore mainly be photospheric molecular ab- 
sorption bands. To search for evidence for the presence of the 
MOLsphere(s) discussed in Sect. l5.ll or for indications of short- 
comings in our model atmospheres, we compare the average nor- 
malizecQresiduals with simulated molecular absorption bands in 
Fig-E] Although our models clearly match the observed spectra 
very well (the better part of the residuals does not reach a 5 % 
level), we detect strong additional absorption by CO and maybe 
also some water and OH. The latter two discrepancies can prob- 
ably be solved by more detailed fine tuning of the atmosphere 
models. Cool CO appears to be very abundant also above the 
photosphere, which is in line with a postulated crucial role in the 
wind-driving mechanism. The fact that no massive MOLspheres 
containing water are observed clearly sets the RSG apart from 
the AGB stars. This finding supports our assumption (Sect. 15. il l 
that the excess at 6 pm and the slope issue at 3.5pm is due to 
some continuum emission source and not an extra-photospheric 
water column. 



6.2. Residual dust features 

The lower panel of Fig.|4]shows the average residuals in the dust- 
dominated part of the spectrum. Many discrepancies are clear: 



8 Normalization of each observed and model spectrum was done by 
inverse scaling with the average flux (F y ) between 3.5 and 4pm. The 
spectra were then subtracted, and the residuals calculated in this way 
averaged over the sample. 
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Fig. 4. Upper panel: The average normalized residuals in the 
blue part of the ISO-SWS spectrum together with simulated 
molecular absorption bands. For more det ails on the computa- 
tion of the molecular templates we refer to lCamil d2002l) . Lower 
panel: The wavelength range of the residuals (black solid line) 
which is dust dominated. 



(1) a strong absorption band starting at 8 pm and possibly ex- 
tending up to 12 pm if we interpret the emission bump at 9.7 pm 
as an underestimated olivine dust fraction, (2) a broad 13 pm fea- 
ture and (3) a broad emission feature at 18 pm. Concerning the 
first problem: a significant SiO layer just above the photosphere, 
which deepens the absorption band, could be present as such 
a layer is d e tected around aOri with MIDI interferometry by 
IPerrin etaT] d2007l) . but the SiO bandhead does not correspond 
very well with the observed residuals around 8 /mi. It may also 
be that the current dust model does not yet contain the right sili- 
cate dust composition/shape/grain size. 
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Source 


Spectral 


T ff 


Melilite 


Olivine 


Alumina 


Men iFenoO 


Carbon 


R 


M 

mod el 


M 


L_/L 




TvDe 

it? 


[K] 


%-mass 


%-mass 


%-mass 


%-mass 


%-mass 


[R*] 


[M yr 1 ] 


[M Q yr -1 ] 




NR VUL 


K3Iab 


4015 


0.56 


0.19 


0.24 


0.00 


0.013 


10 


3.1e-07 


2.7e-07 


0.074 


SIG OPH 


K3Iab 


4015 




















HR 7475 


K4Ib 


3900 




















IRC +40427 


Ml Tab 


3745 


0.43 


0.00 


0.55 


0.00 


0.018 


15 


l.le-07 


5.0e-08 


0.011 


HD 14580 


Ml lab 


3745 




















HD 14404 


Ml lab 


3745 


0.48 


0.00 


0.48 


0.04 


0.007 


15 


1.5e-07 


7.7e-08 


0.012 


ALF SCO 


M1.5Iabb 


3710 


0.74 


0.09 


0.12 


0.00 


0.049 


13 


4.1e-08 


3.3e-08 


0.009 


MU. CEP 


M2Ia 


3660 


0.31 


0.27 


0.41 


0.00 


0.020 


18 


2.9e-07 


7.9e-07 


0.048 


HD 90586 


M2Iab/Tb 


3660 


0.57 


0.14 


0.23 


0.05 


0.011 


10 


1 .4e-07 


8.5e-09 


0.020 


ALF ORI 


M2Iab 


3660 


0.64 


0.16 


0.20 


0.00 


0.001 


13 


6.3e-08 


2.9e-08 


0.022 


HR 1939 


M2Iab 


3660 


0.82 


0.01 


0.16 


0.00 


0.003 


10 


l.le-07 


6.1e-09 


0.006 


YZ PER 


M2Iab 


3660 


0.60 


0.14 


0.26 


0.00 


0.003 


18 


3.5e-07 


2.7e-07 


0.026 


HD 14826 


M2Iab 


3660 


0.42 


0.05 


0.52 


0.00 


0.005 


15 


1.7e-07 


1.0e-07 


0.010 


HD 14242 


M2Iab 


3660 


0.59 


0.00 


0.34 


0.05 


0.017 


15 


1.8e-07 


7.2e-08 


0.023 


V354 CEP 


M2.5Iab 


3615 


0.57 


0.14 


0.28 


0.00 


0.008 


1 1 


5.6e-07 


9.8e-07 


0.087 


BD+35 4077 


M2.5Iab 


3615 


0.78 


0.03 


0.16 


0.02 


0.010 


15 


2.9e-07 


2.0e-07 


0.029 


AD PER 


M2.5Iab 


3615 


0.59 


0.03 


0.34 


0.02 


0.017 


15 


2.4e-07 


1.5e-07 


0.029 


RW CYG 


M3Iab 


3605 


0.53 


0.17 


0.30 


0.00 


0.002 


25 


1.0e-06 


9.4e-07 


0.069 


PZCAS 


M3Iab 


3605 


0.40 


0.30 


0.30 


0.01 


0.002 


40 


1.7e-06 


5.4e-06 


0.109 


SUPER 


M3Iab 


3605 


0.39 


0.19 


0.39 


0.03 


0.008 


15 


2.6e-07 


3.0e-07 


0.027 


ULAC 


M4Iab:e 


3535 


0.25 


0.44 


0.25 


0.01 


0.053 


10 


3.7e-07 


9.0e-07 


0.154 


RHO PER 


M4II 


3535 




















WPER 


M4.5Iab 


3535 


0.09 


0.27 


0.63 


0.00 


0.009 


20 


7.8e-07 


6.4e-07 


0.085 


S PER 


M4.5Iab 


3535 


0.03 


0.27 


0.68 


0.01 


0.022 


15 


1.8e-06 


2.2e-06 


0.302 


RS PER 


M4Iab 


3535 


0.42 


0.21 


0.31 


0.04 


0.017 


13 


3.8e-07 


4.2e-07 


0.063 


ALF HER 


M5Iab 


3450 




















RCEN 


M5IIevar 


3450 


0.33 


0.06 


0.60 


0.00 


0.009 


15 


3.7e-07 


3.3e-07 


0.028 



6.2.1. A 13 jum feature? 

ISoeck etal.l (l2000) find in their analysis that RSG usually do not 
display the 13 /mi feature as seen in AGB stars. We find in gen- 
eral, after model substraction, a feature around 13 //m, but it is 
much broader than the one seen in the lower-mass counterparts. 
Only a few selected targets display the AGB-like 13 /mi feature, 
but these stars (SU Per and S Per) may be different in other as- 
pects as well: S Per is the only target in our sample showing 
strong water bands, and lSpeck et alj d2000t) already suggested it 
to be more similar to lower-mass semi-regulars. The broad ex- 
cess in the residual spectrum can not be isolated in any of the 
spectra prior to model subtraction. Therefore, it can not be pos- 
itively identified as a component, but might be due to some sys- 
tematic problem with our models. The most likely systematic 
effect to cause such an underprediction would be a lack of ab- 
sorptivity in one of the sets of optical properties that we have 
used. 



6.2.2. The 18//m feature 

The 18 //m feature observed in the residuals is most probably the 
feature typical ascribed to the silicates generating the 9.7 /mi fea- 
ture (the olivines). However, starting from laboratory-measured 
optical constants, we find it impossible to reproduce this peak 
position: either it is blue shifted (olivines) or it peaks to the red 
(melilite). The inclusion of either alumina or MgFeO causes the 
model emission band to be centered at the observed wavelength. 
However, this always causes a double peaked band profile which 
is never observed. We note that this discrepancy is not unique to 
the RSG but is also found in the AGB stars with optically thin 



dust shells as studied by iHeras & Honvl d2005). This problem 
may be related to that at 8 /on. 



6.2.3. Detection of PAHs 

As shown in Fig. [5] we clearly detect PAH emission bands at 
6.3, 7.6, 11.3 and 14.2 /mi in 3 stars: AD Per, IRC 40427, and 
ULac. RSPer only shows a significant PAH feature at 7.6 /mi. 
AD Per and IRC 40427 also show a feature at 12.7 /mi. The 
PA Hs in AD Per, RS Per and IR C 40427 were already detected 
by ISvlvester et al.l d 1994 1 1 998b who relate their presence with 
the dissociation of CO molecules by the UV radiation field of 
the chromosphere, as no PAHs are detected in AGB stars. The 
diffuse ISM is known to show PAH emission features but we are 
confident that the bands observed here originate in the CSE of 
the RSG since 1) the observed bands are as strong in the UKIRT 
observations of lSvlvester et al.l(ll994l) as in our ISO-SWS obser- 
vations, in spite of a different aperture, 2) the UKIRT observa- 
tions used a chopping technique, which should remove a large 
fraction of an extended ISM feature and 3) the relative band 
strenghts do not agree with those of the ISM and the 6. 22 /an 
feature is, in some cases, shifted to longer wavelengths. 

Using the dust shell models presented here, we are able to 
isolate the PAH emission features. This allows us to study the 
ban d shapes and stren gths in a manner similar to Peeters (2002) 
and lHonv et al.l d2001l) . Such an analysis can yield important in- 
formation on the structural properties of the PAH family and give 
clues about their formation conditions. This analysis is left to a 
future study. Perusal of Fig. [5] already shows that the PAH fea- 
tures in our sample span an interesting range both in which bands 
are observed, band strength ratios and band-shapes. 
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Fig. 6. Left panel: The correlation between mass-loss rate and the ratio of dust to photospheric luminosity. The slope is that of a 
linear fit, the Kendall values refer to the Kendall rank correlation value and its significance. The Kendall rank correlation r = 

° *n(n— 1) 

where n is the size of the sample, n c is the number of concordant pairs and «,/ is the number of discordant pairs. The significance 
is a value in the interval [0.0, 1.0] where a small value indicates a significant correlation. The relation is clear using both M gr id 
(grey symbols) and the distance-corrected mass-loss rate (black symbols). Middle panel: Mass-loss rate (corrected for luminosity 
and distance) versus stellar luminosity. Right panel: The luminosity- and distance-independent mass-loss indicator versus stellar 
luminosity. Only stars which are certainly RSG were used to measure these correlations. 
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Fig. 5. The 5 to 15 pm spectra of RSG showing PAH emis- 
sion. The grey error bars correspond to 3 cr. Known PAH feature 
wavelengths ar e indicated with a d ash ed line. These w avelengths 
are taken from lHonv et all d2001l) and lPeetersI (12002b . 



7. Analysis 

7.1. Ldust/U as mass-loss rate indicator 

Column 10 presents the ratio between dust and photospheric lu- 
minosity. This ratio can be used as a distance-independent mass- 
loss indicator, which is not trivial, as it must be checked that the 
variation of Ld us t/L* over the sample is due to different mass- 
loss rates and not only a difference in dust composition, e.g. the 
oxides tend to generate much less emission than the magnesium- 
rich silicates, for a similar dust mass. We show our derived re- 
lation between Ld ust /L* and mass-loss rate, based on our mod- 
elling including the differences in dust composition, in the left 
hand panel of Fig. [6] When assuming an identical luminosity for 



all stars, we find an almost one-to-one relation between mass- 
loss rate and dust luminosity fraction. This shows that Ld us t/L* 
is a good mass-loss indicator. When correcting for distance and 
stellar luminosity, the relation is less steep. This is a consequence 
of the correlation between stellar luminosity and mass-loss rate 
discussed in the next section. 

7.2. A correlation between mass-loss rate and stellar 
luminosity 

The middle panel of Fig. [6] shows a strong correlation between 
mass-loss rate and stellar luminosity and therefore also with stel- 
lar mass. The best-fit linear relation is: 

M = io _I6 - 2±L8 L I ' 89±0 ' 36 

As an increase of the mass-loss rate with luminosity, i.e. size of 
the object, is expected just from scaling arguments, we also show 
the dust luminosity fraction as a function of stellar luminosity 
(right hand panel of Fig. [6]). Interestingly, we find evidence for 
a more efficient wind-driving mechanism in the more luminous 
stars. The best-fit linear relation here is: 

T IT ,n-7-5±1.7,1.06±0.30 

We note also that we do not find a significant trend between 
mass-loss rate and effective temperature. 

7.3. Dust composition as a function of mass-loss rate 

A purely empirical indication of the relation between mass-loss 
rate and dust composition can be obtained by simply plotting 
peak position of the 10-pm complex versus dust luminosity. This 
is done in the upper left panel of Fig. [7] With increasing mass- 
loss rate, the peak position shifts from 11 to 10 p.m. For inter- 
mediate values (e.g. BD+35 4077, see top right in Fig.|2]i, it is 
difficult to assign a single peak wavelength. 

We search for correlations between dust composition and 
mass-loss rate based on our modelling and find the following 
trends with increasing mass-loss rate: 
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Fig. 8. Olivine content as a function of temperature and pressure 
at the inner edge of the dust shell. The lines represent regimes of 
constant mass-loss rate, as derived bv lGail & Sedlmavrl (ri986). 



1 . the mass fraction of olivines increases 

2. the melilite content decreases 

3. the MgFeO content may decrease, but the Kendall test yields 
a low significance 

4. the carbon and alumina content don't show a clear trend 

Since the freeze-out scenario actually predicts different com- 
positions for different temperature/pressure combinations in the 
wind, which of course depend on the mass-loss rate, we show the 
olivine content as a function of temperature and pressure at the 
inner edge of the dust shell in Fig. [8] It is clear from this figure 
that indeed, larger fractions of olivine dust are formed at higher 
wind densities as is predicted by the condensation sequence of 
lTielens| (| 19901). and as is seen in AGB stars with low mass-loss 
rates dHeras & Hony 2005; Le bzelter et"ai1l2006l) . 



8. Conclusions 

We have analyzed the ISO-SWS spectra of 27 red supergiant 
stars, out of which 21 show significant mass loss through a 
dusty wind. The determined properties of the winds around RSG 
share some important characteristics with those around AGB 
stars with relatively low mass-loss rates. 1) The dust has a high 
fraction of "simple" dust species like metal-oxides. 2) The winds 
have a relatively low abundance of silicates. 3) The fraction of 
silicates correlates well with the mass-loss rate and/or the den- 
sity and pressure at the base of th e wind, as pr edicted by the 
canonical condensation sequence of lTielensI (11990b . 

However, in some respects, RSG are not just the heavier 
twins of AGB stars: 1) They show molecular bands only of 
di-atomic molecules (not H2O, CO2 or SO2). 2) The general 
slope of the SED from near-IR to mid-IR wavelengths requires 
a source of continuous opacity which, in the case of RSG, could 
be due to amorphous carbon. 3) PAHs are observed in 4 sources, 
suggesting a strong influence by the chromospheric radiation 
field. 

Besides the correlation observed by Josseli n~& Plezl d2007l) 
between photospheric turbulence strength and mass-loss rate, we 
find that also the stellar luminosity strenghtens the wind. 



Although the correlations and detections presented in this pa- 
per are significant (in a rough statistical sense) we are aware of 
possible 1) degeneracies between the fitting parameters, 2) bi- 
ases due to our particular fitting strategy and model assumptions, 
and 3) biases due to our limited set of dust constituents. 

Major progress will be possible in the near future with the 
2nd generation VL TI instrumentation, in particular MATISSE 
(lLopez et a l. 2006) will allow the quasi-instanteneous imaging 
of many RSG from the near to mid-IR. This will allow us to de- 
termine the actual spatial distribution of the different dust com- 
ponents. 
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